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Abstract: 

We study the amplitude distribution of irregular eigenfunctions in systems with mixed classical 
phase space. For an appropriately restricted random wave model a theoretical prediction for 
the amplitude distribution is derived and good agreement with numerical computations for 
the family of limagon billiards is found. The natural extension of our result to more general 
systems, e.g. with a potential, is also discussed. 
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1 Introduction 



The semiclassical behavior of the eigenfunctions of a quantum mechanical system strongly de- 
pends on the ergodic properties of the underlying classical system. The semiclassical eigenfunc- 
tion hypotheses states that the Wigner function of a semiclassical eigenstate is concentrated 
on a region in phase space explored by a typical trajectory of the classical system. In integrable 
systems the phase space is foliated into invariant tori, and the Wigner functions of the quantum 
mechanical eigenfunctions tend to delta functions on these tori in the semiclassical limit || . On 
the other hand, in an ergodic system almost all trajectories cover the energy shell uniformly, 
and hence the Wigner functions of the eigenstates are expected to become a delta function on 
the energy shell. That this actually happens for an ergodic system for almost all eigenstates 
follows from the quantum ergodicity theorem see [§-0], and 0,|8[] for billiards (the relation of 
the quantum ergodicity theorem to the semiclassical behaviour of Wigner functions is explicitly 
derived for Hamiltonian systems in ||). However, a generic system is neither integrable nor 
ergodic fllPl , but has a mixed phase space in which regular regions (e.g. islands around stable 
periodic orbits) and stochastic regions coexist. Whether these numerically observed stochastic 
regions are ergodic and of positive measure is an open question, see |lT| for a review on the 
coexistence problem. The eigenfunctions in mixed systems are expected to be separated into 
regular and irregular eigenfunctions according to an early conjecture by Percival Jl2| which has 
been numerically confirmed for several systems, see e.g. [P^T5|. In addition, at finite ener- 
gies there is a small (semiclassically vanishing) fraction of "hierarchical states" which are of 
intermediate nature, and localize in regions bounded by cantori Jl7 . 



Beside the localization properties of the Wigner function, also the local amplitude fluctu- 
ations of the eigenfunctions strongly depend on the classical system, as has been pointed out 
in ]TJ]. The basic idea is that an eigenfunction can be represented locally as a superposi- 
tion of de Broglie waves with wavelength determined by the energy, and momenta distributed 
according to the semiclassical limit of the Wigner function. In a chaotic system one therefore 
expects an isotropic distribution of the momenta. If one additionally assumes that the phases 
are randomly distributed, one obtains locally a Gaussian amplitude distribution of a typical 
eigenfunction in a quantum mechanical system with chaotic classical limit. For instance in a 
chaotic billiard a Gaussian amplitude distribution is expected, and this has been confirmed by 
several numerical studies, see e.g. [|i9|-f2~4]l . Furthermore, predictions of the random wave model 
for the maxima of chaotic eigenfunctions have been derived and numerically tested for several 



systems in |Pq| . In contrast, in an integrable system the localization of the Wigner function on 
the invariant tori enforces a more coherent superposition of the de Broglie waves, leading to a 
regular structure of the eigenfunction ]T|| . 

Our aim is to determine the amplitude distribution for irregular states in systems with a 
mixed classical dynamics. We assume that the motion on a stochastic region D in phase space 
is ergodic and that the statistical properties of eigenfunctions can be described by a random 
wave model restricted to D see the following section for a precise definition. The derivation 
shows that locally the fluctuations are Gaussian with a position dependent variance which is 
given by the classical probability density on position space defined by the ergodic density on 
D. Thus the resulting amplitude distribution may be significantly different from a Gaussian. 
In section |3| we compare the theoretical prediction of the restricted random wave model with 
numerical computations. 
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2 Amplitude distribution for the restricted random wave 
model 



In this section we consider a restricted random wave model for two-dimensional Euclidean 
quantum billiards in order to describe the statistical properties of irregular eigenf unctions in 
systems with a mixed classical phase space. The quantum mechanical system is defined by the 
Euclidean Laplacian on a compact domain Q C M 2 with suitable boundary conditions on the 
boundary dQ. (Usually one chooses Dirichlet conditions.) The quantum mechanical eigenvalue 
problem is given by 

-Ai[) n (q) = E n i/; n (q) , with ^ n (q) = for q G dQ , (1) 

and we are interested in the behavior of the eigenf unctions ip n in the semiclassical limit E n — > oo. 

The corresponding classical system is given by a free particle moving along straight lines 
inside the billiard, making elastic reflections on the billiard boundary dQ. The phase space is 
T*Q = IR 2 x Q, and the Hamiltonian is H(p, q) = |p| 2 . Since the Hamiltonian is scaling we can 
restrict our attention to the equi-energyshell with energy E = 1, 

S*Vl :={(p,q)GR 2 xfl; |p| = 1} . (2) 

Introducing polar coordinates (r, </>) for the momentum p, we can parametrize S*Q by (0, q) G 
[0, 2tt) x Q where <ft is the direction of the momentum. In these coordinates the Liouville measure 
on S*Q is given by 

d/i = d0d 2 g, (3) 

which is invariant under the Hamiltonian flow on S*Q. 

Now let D C S*Q be an open domain which is invariant under the classical flow, and on 
which the flow is chaotic. The existence of such a domain where the flow is, for instance, 
ergodic, is an open problem. But numerically one observes invariant domains on which the 
flow is at least irregular in the sense that most orbits are unstable, and regular islands inside 
this domain are very small. The uncertainty principle implies a finite quantum mechanical 
resolution of phase space quantities at finite energies. Therefore at finite energies the small 
islands of such an irregular domain are not resolved by the quantum system. 

So we expect, in the spirit of |], that the statistical properties of irregular eigenfunctions 
associated with D can be described by those of a superposition of plane waves with wave 
vectors of the same lengths and directions distributed uniformly on D. If we furthermore 
assume random phases we arrive for real valued functions at the following restricted random 
wave model, which is a superposition of plane waves of the form 



4:71 



-0RRWM,D(q) = J vo ^ D ^ N Xd(K, q) cos(k n q + e r 



(4) 



Here Xd(") * s the characteristic function of D, the phases e n are independent random variables 
equidistributed on [0, 2tt], and the momenta k n G M 2 are independent random variables which 
are equidistributed on the circle of radius \[E. So the characteristic function Xd(-) ensures the 
localization on D. Furthermore, it is natural to take N ~ \Z~E~, the scaling of the number of 
line segments of a typical Heisenberg-length orbit. By vol(D) the volume of D measured with 



3 



the Liouville measure (^) is denoted. With this choice of normalization the expectation value 
of the norm | |^rrwm,d| | is one. 

Let us first consider the value distribution -P q (V ; ) of V'RRWM.i^q) at a given point qGfi. Our 
restricted random wave model (f|) is a sum of identical independent random variables which 
have zero mean and whose variance is given by 



a 2 (q) = E ( vo jp) (X£>(kn»q) cos(k„q + e n )) 2 ^ 

r 

Xr»(e(0),q) d</> , 



2vr 

1 



vo\(D) 



(5) 



where e(0) := (cos(0), sin(0)) denotes the unit vector in ^-direction. So by the central limit 
theorem we obtain for E — > oo, i.e. N — > oo, a Gaussian distribution of i/j-rjwm,d(<1) at q, 



with variance given by (||). If the classical dynamics on D is ergodic, then the variance cr 2 (q) 
is exactly the probability density of finding the particle in the point q G Q if it moves on a 
generic trajectory in D. So cr 2 (q) is the classical probability density in position space. 

By integrating eq. (^) over Q we obtain the complete amplitude distribution as a mean over 
a family of Gaussians with variances given by (|5|), 

Prrwm,dW = [ P M d 2 q (7) 



vol(fi) 




vol(fi) 7 V 2vra 2 (q) ^ V 2(X 2 (q) 



exp -t^--^ 2 d 2 g . (8) 



So the amplitude distribution is completely determined by the classical probabality density (0), 
and it will be typically non Gaussian if cr 2 (q) is not constant. 

The moments of the distribution (||) can be computed directly and turn out to be propor- 
tional to the moments of the classical probability density <r 2 (q), 

J ^Prrwm,!) W d^ = P 2k J [a 2 (q)] fc dg , (9) 

n 

where the factor p 2 fc = denotes the 2/c'th moment of a Gaussian. The odd moments are of 
course zero. Note that the second moment is always l/vol(fl), due to the normalization of ip. 

If the system is ergodic one has cr 2 (q) = vol ^ and we get the classical result that -Prrwm,d(' ? /') 
is Gaussian with variance a 2 = vol "^ • However, if cx 2 (q) depends on q then the corresponding 
distribution can show deviations from the Gaussian distribution. In particular, if <r 2 (q) = 
for some region fl' C fi, we get a contribution ^^ (^) to the corresponding distribution of 
-Prrwm,d(V ; ) as the integrand in (|7|) tends to a 5 distribution as cr 2 (q) — > 0. 

Finally, we would like to point out that the main ingredient in the formula (|7|) is the as- 
sumption that the local amplitude distribution of an irregular eigenfunction around a point q in 
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position space is Gaussian, with a variance given by the classical probability density in position 
space cr 2 (q), defined by the projection of the invariant measure on D to position space. Clearly 
this assumption is not restricted to billiards, but is expected to be true for arbitrary quantum 
mechanical systems for which the underlying classical system contains chaotic components in 
phase space. So the formula (0) is expected to be valid in far more general situations, with 
<r 2 (q) denoting the classical probability density defined by the ergodic measure on the chaotic 
component. 



3 Comparison with irregular eigenfunctions 

We now compare the predictions of the restricted random wave model with the results for 
some numerically computed eigenfunctions. As systems to study the amplitude distribution of 
irregular states in mixed systems we have chosen the family of limagon billiards introduced by 
Robnik P6|,p7| with boundary given in polar coordinates by p{tp) = 1 + £:cos(y?), <p G [— 
with e G [0, 1] being the system parameter. We consider the case e = 0.3, for which the 
billiard has a phase space of mixed type [2B|, see fig. [TJ. In E8[] examples of eigenstates far 
in the semiclassical regime have been studied in this system, and in particular the amplitude 
distribution has been studied numerically, but no analytical predictions have been made. 

First we have to determine the classical position space probability density a 2 (q) of the 




Figure 1 : Plot of several stable and irregular orbits in the Poincare section V of the limagon billiard 
for e = 0.3. Here V is parametrized by the (rescaled) arclength coordinate s G [—4, 4] along the 
billiard boundary and p £ [—1,1] which is the projection of the unit velocity vector on the tangent 
in the point s after the reflection. 
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ergodic measure on the invariant domain D. The normalized ergodic measure on D is given by 



dfi D {<f),q) = .) n ^ o(e^),q) d<pd 2 q , 
vol{D) 

so we can express the variance cr 2 (q) as a mean value 

^ 2 (q)= f a(q-q') d^(0',q') . (10) 



S*f2 

As the motion on £) is assumed to be ergodic, in order to determine cx 2 (q) we could replace 
the integral over S*Q by a time average over a typical trajectory of D and the 5 function 
by a smoothened 5 function, e.g. a narrow Gaussian. However, as we will see below, the 
eigenfunctions turn out to be not concentrated on the whole chaotic component, but rather on 
a subset which is almost invariant in the sense that it is bounded by partial barriers in phase 
space. Since at finite energies quantum mechanics has only a finite resolution in phase space, 
these partial barriers appear like real barriers. But since any classical trajectory will pass such 
a barrier after a certain time, the time average is not suitable for the determination of cr 2 (q) 
in such a situation. 

For a more direct approach to determine <r 2 (q) we use the Poincare section V = {(s,p) ; s G 
[—4,4], p G [—1,1]}, which is parametrized by the (rescaled) arclength coordinate s (corre- 
sponding to <p G [— 7T, 7t]) along the boundary dQ and the projection p of the unit velocity 
vector on the tangent in the point s after the reflection. Let T> C V be the projection of the 
region D in the energy shell S*Q := {(p, q) G 1R 2 x Q ; ||p|| = 1} to the Poincare section. This 
projection is defined as follows: to a point (e(0), q) G D we can associate the trajectory which 
passes through q in direction e(0), then s(<f>, q) is defined as the first intersection with the 
boundary dQ when traversing the trajectory backwards from q, and p(<j>, q) := e(0)T(s(0, q)) 
which is the projection of the unit velocity vector e(0) to the unit tangent vector T(s(0, q)) to 
dQ at s(<p, q). 

For a given point q we therefore get a curve parameterized by <fi 

(p(0,q),s(0,q))GP. (11) 
Since x^(e(0),q) = Xv((p(<(>, q), s(0, q)), we get 

^ 2 (q) = ^j|^y y X^((p(0,q),s(0,q)) d0 , (12) 
o 

and therefore we have to determine the fraction of the angular interval(s) for which the curve 
(|i~ID is in T>. That is, one has to determine the angles 0i ntry (q) and 0i Xlt (q) where the curve 
(jTTT) enters or leaves the region "D, i.e. the intersection points of (|i"T|) with the boundary of T>. 
In terms of these angles we obtain 

^ 2 (q) = ^rj) E ^ xit (q) - C try (q) , (is) 

which is proportional to the fraction of directions in the ergodic component visible from the 
point q. 
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With this classical probability density cr 2 (q) one can compute the corresponding amplitude 
distribution via eq. ||. If <r 2 (q) = for some region, then the local amplitude distribution @ 
becomes a delta function, and it is necessary to consider for a concrete comparison a binned 
distribution, 



ip+Aip/2 



Alp 



P(ip') dip' 



ip-Atp/2 



2IQI 



erf 



Mq) J 



erf 



ip - Aip/2 



d 2 g . 



(14) 
(15) 



We now use a Husimi Poincare section representation of the eigenstate (see e.g. p9|,|3T 
to determine the boundary of the relevant component I? by a spline approximation. The 
Poincare Husimi representation of an eigenfunction ip n in a billiard is defined by projecting the 
normal derivative u n (s) of an eigenfunction VVi(q) at the boundary onto a coherent state on 
the boundary. The coherent states, semiclassically centered in (s,p) G V, are defined as 



\ — J exp (ipk (s' — mL — s)) exp ( — — (s' — mL — s) j 

\ / <m — — rv^ ^ * 



(16) 



where s' e [—4, 4] , a > and L = 8 is the total (rescaled) length of the boundary. This 
definition is just a periodized version of the standard coherent states, 
function of a state ip n with normal derivative u n (s) is then defined as 



The Poincare Husimi 



H n {s,p) 



2tt 



\u. 



ds 



J (s,p),k n 



U,. 



s') ds' 



(17) 



with k n = \[E^\ the prefactor ensures the normalization J J H n (s,p) dpds = 1. 

A first example is shown in fig. [|. In a) a high lying eigenfunction (E = 1002754.70 . . . , 
approximately the 130568 th state of odd symmetry) in the limagon billiard with e = 0.3 is 
shown as density plot (black corresponding to high intensity of \ip\ 2 )- In b) the corresponding 
Husimi representation on the Poincare section is shown. The boundary of the irregular region 
T> is described by a cubic spline which is shown as full curve. With these boundary curves we 
can use (|HJ) to compute cr 2 (q), which is shown in fig. ^|c). Finally, in fig. |2]d) the comparison 
of the amplitude distribution of ip with the prediction of the restricted random wave model 
is given. Clearly, P(ip) is non-Gaussian, and the agreement is very good. Table [I] lists the 
first moments and also very good agreement of the results using (^[) and the moments of ip is 
found. Both the resulting amplitude distribution Prrwm,d and the moments turn out to be 
quite robust with respect to small changes of the selection of T>. Notice that we have rescaled 
cr 2 (q) such that the variance of the distributions is one. 

Another example is shown in fig. The eigenfunction (E = 1003030.75 . . . , approximately 
the 130607 th state of odd symmetry) plotted in a) has a quite large region in the center where 
it is almost vanishing. So from this alone the amplitude distribution is expected to show a very 
clear deviation from the normal distribution. Using the same procedure as in the previous case 
we determine V, compute cr 2 (q) and then Prrwm(^)- The comparison of the prediction with 
P(ip) is shown in fig. |d). The strongest deviation occurs for ip w 0. The peak of Prjwm(iP) 
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Figure 2: In a) a high lying eigenfunction (E = 1002754.70 . . . , approximately the 130568 state) 
in the limagon is shown as grey scale plot (black corresponding to high intensity). In b) the 
corresponding Husimi function on the Poincare section is shown together with the boundary (full 
curves) of the region on which the eigenfunction is concentrated. In c) a density plot of 0" 2 (q), 
computed via eq. (|l^), is shown. In d) the cumulative amplitude distribution of the eigenfunction 
is compared with the prediction of the restricted random wave model (RRWM); on this scale no 
differences are visible. The left inset shows P(ip) and for the right inset a logarithmic vertical scale 
is used to emphasize the tails of the distribution. For comparison the normal distribution is shown 
as grey curve. g 
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Figure 3: The same plots as in the previous figure are shown for another high lying eigenfunction 
(E = 1003030.75. . . , approximately the 130607 th state). In this case there is a deviation of the 
amplitude distribution of the eigenfunction from the prediction of the restricted random wave model 
around ip = 0. This is because <r 2 (q) = in the central region, whereas the eigenfunction does not 
vanish there (see the text for further discussion) . For the tails of the distribution the agreement of 
the two distributions is again very good. 



9 



IIlUIIlcIll 


Example 
eigenfunction 


1, fig- I 
RRWM 


Example 
eigenfunction 


2, fig- 1 
RRWM 


T"l f"\Y*TYl Q 1 riict"ri r\n"f"1/~\T"l 


4 


4.39 


4.46 


3.85 


3.75 


3 


6 
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269 


269 
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10 


2199 


2501 


3774 


3841 


945 



Table 1: Comparison of the even moments for the distributions of the eigenfunction and the 
restricted random wave model (RRWM) eq. (||) . The last column lists for comparison the moments 
of the normal distribution. 
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Figure 4: For the three domains indicated in the inset the local amplitude distribution is shown 
(for the same state as in fig. |||). We observe very good agreement with the expected Gaussian 
behaviour (shown as dotted curves) with position dependent variance. Notice that the non-zero 
width of the distribution for the region B corresponds to the widening of the ^-contribution, see 
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at ip — is due to the fact that <r 2 (q) = for the region in the center of the billiard. The 
eigenfunction, however, is not exactly zero, but shows a decay into that region and thus still 
fluctuates there. This causes a broadening of the 5-contribution, which is clearly visible in the 
plot of P(ip) in fig. |3]d). For > 0.25 this region is not relevant anymore, and the agreement 
of P(ip) and -PrrwmC^) is ver Y good. In the right inset of fig. |]d) the distribution is shown with 
a logarithmic vertical scale to illustrate the agreement of the distributions even in the tails. 

The moments, computed via eq. @, are listed in table |l|. The agreement of the moments 
of the eigenfunction with the prediction of the restricted random wave model is quite good. 
All moments of the two examples are larger than those of a Gaussian, corresponding to the 
larger tails. Compared to the moments of the restricted random wave model those of the 
eigenfunctions tend to be smaller, in particular for the larger moments. This is reasonable, 
as an actual eigenfunction is always bounded, which reduces higher moments compared to the 
result of eq. @. 

We furthermore have tested our basic assumption @, that the local value distribution of 
a sufficiently high lying eigenfunction is Gaussian with a variance given by the local classical 
probability density associated with D, more directly. To this end we have computed the value 
distribution of the eigenfunction in fig. ^| for three small regions on which <r 2 (q) is almost 
invariant, and we therefore expect a Gaussian. The results are shown in fig. |], and good 
agreement with the prediction (^) is found. Since many fewer wavelengths are contained in 
these small domains than in Q the statistics is of course not as good as for the full system, 
but the results give strong support for a local Gaussian behavior. The variances for the two 
domains A and C coincide with the expected classical one cr 2 (q). But for domain B the observed 
variance is larger than cx 2 (q) = 0. This corresponds to the widening of the delta peak in fig. |||, 
and is due to the fact that the eigenfunction cannot become exactly zero on some open set at 
finite energies, but instead fluctuates around zero. 



4 Summary 

In this work we have extended the random wave model for eigenfunctions from the case of 
chaotic systems to the case of irregular eigenfunctions in systems with mixed phase space. Our 
main result is one particular prediction of this model, namely the amplitude distribution (]?]) 
of irregular eigenfunctions. Numerical tests have been performed for two high lying eigenfunc- 
tions of the limagon billiard with e = 0.3, and impressive agreement, even in the tails of the 
distribution, with the theoretical prediction was found. 

The physical picture underlying our analysis is that the local hyperbolicity in the irregular 
part of the phase space forces the eigenfunctions localizing on this part of phase space to behave 
locally like a Gaussian random function with a variance given by the classical probability density 
in position space defined by the uniform measure on the irregular component. Taking the mean 
over all these local Gaussians with varying variance gives our result for the global amplitude 
distribution. We have tested this intuitive picture by computing local amplitude distributions. 
The agreement of these with the Gaussian prediction is very good, giving further strong support 
to the picture of local Gaussian fluctuations with variance determined by the underlying classical 
system. A natural further question relates to the correlations of such eigenfunctions between 



different points in position space; this topic is addressed in |31 . 

Although we have restricted our study to Euclidean billiards, the general picture of local 
Gaussian fluctuations is of course not limited to these special type of systems. We therefore 
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expect our results to be valid for irregular eigenfunctions in arbitrary systems (e.g. systems 
with potential), with o" 2 (q) defined as the projection of the ergodic measure on the irregular 
component to the position space. 
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